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Abstract 

We review and assess a part of the recent work on Casimir apparatuses in the weak gravitational 
field of the Earth. For a free, real massless scalar field subject to Dirichlet or Neumann boundary 
conditions on the parallel plates, the resulting regularized and renormalized energy-momentum 
tensor is covariantly conserved, while the trace anomaly vanishes if the massless field is conformally 
coupled to gravity. Conformal coupling also ensures a finite Casimir energy and finite values of the 
pressure upon parallel plates. These results have been extended to an electromagnetic field subject 
to perfect conductor (hence idealized) boundary conditions on parallel plates, by various authors. 
The regularized and renormalized energy-momentum tensor has been evaluated up to second order 
in the gravity acceleration. In both the scalar and the electromagnetic case, studied to first order 
in the gravity acceleration, the theory predicts a tiny force in the upwards direction acting on 
the apparatus. This effect is conceptually very interesting, since it means that Casimir energy is 
indeed expected to gravitate, although the magnitude of the expected force makes it necessary to 
overcome very severe signal-modulation problems. 
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I. INTRODUCTION 



Ever since Casimir discovered that suitable differences of zero-point energies of the quan 



several 



tized electromagnetic field can be made finite and provide measurable effects 
efforts have been produced to understand the physical imphcations and applications of this 
property . In particular, we are here going to review the recent theoretical discovery 

that Casimir energy gravitates In Ref. 9], this was proved as part of an investi- 

gation that led, for the first time, to the evaluation of the energy-momentum tensor of a 



Casimir apparatus in a weak gravitational field (cf. the work in Ref. Uj). In that piece of 
work. Maxwell theory was quantized via functional integral, with perfect conductor bound- 
ary conditions on parallel plates at distance a from each other. On using Fermi-Walker 
coordinates, where the {xi,X2) coordinates span the plates, while the z = X3 axis coincides 
with the vertical upward direction (so that the plates have equations z = and z = a, re- 
spectively), and working to first order in the constant gravity acceleration g, the spacetime 
metric reads as jsl 

ds"^ = -c^ {l + £^ dt^ + dxl + dA + dz^ + Oda;^, (1.1) 

where e = 

Our paper provides a review of some key findings by the authors and by other research 
groups interested in the same topics. For this purpose. Sec. II studies the Feynman Green 
function for the scalar wave operator to zeroth and first order in e. Sec. Ill obtains the 
resulting regularized and renormalized energy-momentum tensor while Sec. IV evaluates 
Casimir energy and pressure upon the plates. All of this with Dirichlet conditions on the 
plates for the Green function. The case of Neumann boundary conditions is considered in 
Sec. V, while the electromagnetic analysis is summarized in Sec. VI. Concluding remarks 
are presented in Sec. VII. 

II. FEYNMAN GREEN FUNCTION TO ZEROTH AND FIRST ORDER 

To first order in the e parameter of Sec. I, the only nonvanishing Christoffel symbols 
associated with the metric (1.1) are 
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We now compute the wave operator the Feynman Green function of the hyperbohc 
operator — C,R), and eventually the Hadamard function and the regularized energy- 
momentum tensor. 

Indeed, a Green function of the scalar wave operator obeys the differential equation 

{n-mG{x,x') = ~^-^. (2.2) 

The Feynman Green function Gi;' is the unique symmetric complex- valued Green function 
which obeys the relation 

6G = G6F G, 

where F is the invertible operator obtained from variation of the action functional with 
respect to the field. This definition is well suited for the purpose of defining the Feynman 



Green function even when asymptotic flatness does not necessarily hold [12|. 

In our first-order expansion in the e parameter, the scalar curvature gives vanishing 
contribution to Eq. (2.2), which therefore takes the form (hereafter Q° = rj^^'^d^du) 

(□o + _:L^^ + r'„„^^|-)G(.,x') = -^. (2.3) 

V {a + ez)dt^ ^"{a + ez)dzj ' 



We now follow our work in Ref. |9| and assume that the Feynman Green function admits 
the asymptotic expansion 

Gf[x, x') ~ G(°)(x, x') + £GW(x, x') + 0(£2). (2.4) 

Its existence is proved by the calculations described hereafter. Indeed, by insertion of (2.4) 
into (2.3) we therefore obtain, picking out terms of zeroth and first order in e, the pair of 
differential equations 

□ °G'(°)(a;,a;') = /°Ha;,a;'), (2-5) 

^''G^^\x,xJ) = J^'-Xx.x'), (2.6) 

having set 

j(°)(x,a;') = -5(x,x'), (2.7) 

,a.,,,,.)^±„,,,0-(£^ + i-|)G<«'(.,x'). (2.8) 

Our boundary conditions are Dirichlet in the spatial variable z. Since the full Feynman 
function Gp{x^x'^ is required to vanish at 2; = 0,a, this implies the following homogeneous 
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Dirichlet conditions on the zeroth and first-order terms: 



= 0, (2.9) 

2=0, a 



G^^\x,x')\ =0. (2.10) 

I 2=0, a 

To solve Eqs. (2.5) and (2.6), we perform a Fourier analysis of G^^^ and G^^\ which remains 
meaningful in a weak gravitational field j^, by virtue of translation invariance. In such an 
analysis we separate the z variable, i.e. we write (cf. joj]) 

GW(x,x') = [ ^!^!^^(o)(^,^')e^fc^^-(-^x-xY)-i/cO(.o-xi)^ (2.11) 
J (27r)'^ 

and similarly for G^^\x,x'), with a "reduced Green function" 'y^-^\z,z') in the integrand 
as a counterpart of the zeroth-order Green function 7^(2, z') in (2.9). Equations (2.3) 
and (2.4) lead therefore to the following equations for reduced Green functions (hereafter 

+\A^^'\z,z') = -Siz,z'), (2.12) 



dz'^ 



By virtue of the Dirichlet conditions (2.9), 7*-°'' reads as 

A sin(Aa) 

where ^< = mm(z,z'), z> = max(z, z'). The evaluation of the reduced Green function 7*^^) 
is slightly more involved. For this purpose, we distinguish the cases z < z' and z > z', and 
find the two equations 

(^-^ + X'y^\z,z')=£\z,z'), (2.15) 

where 

.(1) 1 Acos(Az) -2z/cgsin(A2;) . -r ^ / fo Ta\ 

J- =7, X ■ fx \ sm(A(2; - a)) if 2; < z , (2.16) 

2a Asm(Aaj 

^ _^Acos(A(z-a))-2zfcgsin(A(z-a)) ^ 

2a Asin(Aa) ^ ^ ^ ^ 

We have therefore two different solutions in the intervals z < z' and z > z' . In this case the 
differential equation (2.15) is solved by imposing the matching condition 

^^yiz',z')=^^^\z',z') (2.18) 



jointly with the jump condition 



dz 



7+ 



_ ^^(1) 



z 

2^' 



(2.19) 



z' 9z 

Equation (2.18) is just the continuity requirement of the reduced Green function 'j^^\z,z') 
a.t z = z', while Eq. (2.19) can be obtained by integrating Eq. (2.13) in a neighborhood of 



z , smce 



lim -^7*-^'' 

e-O dz ' 



z z' 

lim / — 6(z,z')dz = — . 

2a ^ ^ 2a 



(2.20) 



Bearing in mind Eq. (2.14) we can therefore write, for all z, z' , 



7«(^,^') 



4aA2 



ikl-X^)iz + z')-kl(z^l^+z'^-^ 



l^'\z,z') 



,9 aSin(Az) sin(Az') 



sin (Aa) 



(2.21) 



III. REGULARIZED AND RENORMALIZED ENERGY-MOMENTUM TENSOR 

In the previous section we have focused on the Feynman Green function Gp because it is 
then possible to develop a recursive scheme for the evaluation of its asymptotic expansion 
at small e. However, we eventually need the Hadamard function H{x,x'), which is obtained 



as 
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H{x,x') = 2ImGF(a;,x') ~ 2lm{G'^^\x,x') + eG^^\x,x')) + 0{e^). 



(3.1) 



The coincidence limits in the formula of the regularized and renormalized energy-momentum 
tensor make it necessary to perform the replacements 

where P^, is the parallel displacement bivector [l^ 

PC, ~ diag (l + ^(^' - z), 1, 1, l) + 0{e'). (3.3) 

Hence we get the asymptotic expansion at small e of the regularized energy-momentum 
tensor according to (hereafter we evaluate its covariant, rather than contravariant, form) 



(r,.)~(r(°)) + e(T«) + 0(5^) 



(3.4) 
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where, on defining s = nz/a, s' = irz'/a, we find 



and 



IT 



— lim 



TT 



/l o\ 

0-100 
0-10 



3 — 2 sin s 

'• 4 

sm s 













3/ 










o\ 







-1 
















-1 




















1440a4 sin'^ s 



311 637 1 „^ , ^ 

TT s H (437r - 81s) cos 2s 

40 40 10^ ^ 



s — 37r 

+ — — — cos 4s + 5 sin 2s + 2{tt — s)s(sin 2s — 6 cot s) 



40 



+ U-7 



vr 



6 J 48a4 sin^ s 
+ (tt — s)s(sin2s — 6 cot s 



2(7r + s)(2 + cos 2s) + - sin 2s 



lim 



TTS 



TT 



TT - 2S + 



sin^ s 



(2(^ - 2s) (-2 + 
V \ sm s / 



7200a4 

+ cot s ^5 + 2(7r — s)s — 6(7r — s) . ^ 



TT 



sm S/ 
(ll(7r — s)s — 1) cos s 



6 y 960"^ sin^ s 
+ ((tt — s)s + 1) cos3s — 2(7r — 2s)(3sins + sin 3s) 



TTS / 1 



TT 



TT COS S 



(3.5) 



(3.6) 



(3.7) 

(3.8) 
(3.9) 



1440a4 ' 720a4 ' V"" 6 J 160^ sinks' 
The next step of our analysis is the renormalization of the regularized energy-momentum 
tensor. For this purpose, following our work in Ref . , we subtract the energy-momentum 
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tensor evaluated in the absence of bounding plates, i.e. 

/l o\ 
0-100 



lim - , , 



and 



lim 



TTS 



0-10 
\0 3/ 



/l o\ 





\0 0/ 



(3.10) 



(3.11) 



To test consistency of our results we should now check whether our regularized and renormal- 
ized energy-momentum tensor is covariantly conserved, since otherwise we would be outside 
the realm of quantum field theory in curved spacetime, which would be unacceptable. In- 
deed, the condition 

V^{T,.) = (3.12) 
yields, working up to first order in e, the pair of equations 

d 



dz 



{ni^) = Q, (e^term) 



|(T«) + ^(rf)) + (Ti?))=0(.Uerm), 

which are found to hold identically for all values of ^ in our problem. 

The trace of {T^^) is obtained as 



(3.13) 
(3.14) 



r 



(3.15) 



from which we find a .^-dependent part 

37r2(2 + cos2s) 



TC 



■^a^ sin^ s 



32a^ sin^ s 

(1 + (tt — s)s) COS 3s + 2{n — 2s)(3sins -|- sin 3s 



1 — ll(7r — s)s) cos s 



(3.16) 



Interestingly, the value ^ = | which yields conformal invariance of the classical action is the 
same as the value of ^ yielding no trace anomaly [3] • 
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IV. CASIMIR ENERGY AND PRESSURE 



In order to evaluate the energy density p of our "scalar" Casimir apparatus, we project 
the regularized and renormalized energy-momentum tensor along a unit timelike vector 



0,0,0 . This yields 



IT 



+ 



IT 



1440a4 7200a4 



Svr + 6s 



10 



sin^ s 



2(71 - 2s) 



X 



3 \ / , , , S(— TT + S] 

-2 + — ^ I + cot s (5 + 2(7r - s)s + 6- 



+ U 



sm^ s 



1\ r7r2(2 + cos 2s) 
6/1 8a^ sin^ s 



sin^ s 



TC 



192a4 sin^ s 



5 + 22(7r — s)s ) cos s 



+ ^5 + 2(7r — s)sj cos 3s — 4(7r — 2s) (3 sins + sin 3s) e^. 
The energy E stored within our Casimir cavity is given by 



E 



-gp, 



(4.1) 



(4.2) 



where d^T, is the volume element of an observer with four-velocity u^, and Vc is the volume 
of the cavity. The integration used here requires the use of approximating domains, i.e. the 
^-integration is performed in the interval {C,a — (), corresponding to ^{(,a — () in the s 
variable, taking eventually the C ^ limit. We thus obtain 131] 



E, 



1440a3 5760a3 



6 / Aa? 



4/ C^o sin'^C' 



(4.3) 



where A is the area of parallel plates. Note that the conformal coupling value ^ 



IS 



picked out as the only value of ^ for which the Casimir energy remains finite. In this case, 
reintroducing the constants c and writing explicitly £, we find {l3] 



hc-n^ A ( Iga 



1440 a3 

In the same way, the pressure on the parallel plates is found to be [1 

/ l\ -ne coss 

PJz = 0) = J H T - f - - J lim — 5-, 

' 480a4 1440a4 V 6 / 16a4 s^o sin^ s 

2 



vr 



+ 



480a4 1440a4 



Tie , . cos s 
lim 



6 / 16a^ s^TT sin s 



Once again, one can get rid of divergent terms by setting ^ = |, which leads to 



Pc{z = 0) 



tt"^ he / 2 ga 
480^ V^3^ 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



Pc{z = a) 



vr^ he / 2 ga 
'480 V ^3c2" 



(4.8) 



To obtain the resulting force one has to multiply each of these pressures by the redshift r of 
the point where they act, relative to the point where they are added i.e., 



|fi'oo(-Pa( 



1 + ~ ^q)' 



I fi'Ool -Padded) 

to leading order in ^. Thus, a net force is obtained of magnitude 



(4.9) 
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480c2 



Z2 — Zi 



^9 



1440c2 



Z2 — Z\ 



TT^ Ahg _ \El\ 
1440^ ~ 



(4.10) 



having defined E'^ = —j^hc-^, which points upwards along the 2;-axis and is in full agree- 
ment with the equivalence principle. 

V. NEUMANN BOUNDARY CONDITIONS 



When the reduced Green functions obey instead Neumann boundary conditions on par- 
allel plates, i.e. 

' - ^ -0, V2 = 0,l, (5.1) 



dz 



^7 



z=Q 



dz 



our work in Ref. \)JS\ has found, by an analogous procedure, the regularized and renormalized 
energy-momentum tensor to first order in e, with trace 



r^^9"'{T,u) 

TT 1 



6/ 32a4sin^s 



67r(3 sin s + sin 3s) 



— e 



1 + — s)s) coss — (1 — (tt — s)s) cos 3s 
2(7r — 2s)(3sins + sin3s) |, 



(5.2) 



which vanishes in the case of conformal coupling, as with Dirichlet boundary conditions. 
Moreover, the Casimir energy stored between the plates, the pressure on parallel plates 
and the net force are then found to agree completely with (4.4), (4.5)-(4.6) and (4.10), 
respectively. 
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VI. ELECTROMAGNETIC FIELD 



The work in Ref. [lb] has instead exploited the fact that, to first order in the small 
quantity gz, the line element (1.1) coincides with the Rindler metric 

-] dt^ + d^'^ + dxl, (6.1) 

where ^ = ^ + z = $,i+z. In a fully covariant analysis of Feynman Green functions in Rindler 
spacetime, the components of the Maxwell energy-momentum tensor have been evaluated 
up to second order in g with perfect conductor boundary conditions, finding that, as 2; — > 0, 

(0|T/|0) = (|T/|0)~-^ + O(.-). (6.4) 

Since T^z is now found to diverge on approaching the plates, no definite meaning can be 
given to the gravitational correction to the Casimir pressure. Moreover, the divergences in 
T^* are such that the resulting correction to the total mass energy of the cavity is infinite, 
even on taking the principal- value integral of T/. 



VII. CONCLUDING REMARKS 

The literature on the behaviour of rigid Casimir cavities in a weak gravitational field 
predicts, on theoretical ground, that Casimir energy obeys exactly the equivalence principle, 
and hence the Casimir apparatus should experience a tiny push in the upward direction. The 



formula for t 
over modes 
9 



le push has been obtained in three different ways, i.e. a heuristic summation 



16 



171, or a variational approach 

n 



7 



, |8l], or an energy-momentum analysis 



15| . Moreover, the work in Ref. [18|] has shown that Casimir energy for a configuration 
of parallel plates gravitates according to the equivalence principle both for the finite and 
divergent parts. This suggests that such divergent parts can be absorbed by a process of 



renormalization 



18|. 



It now remains to be seen whether this interpretation is viable in all configurations of 
physical interest. Moreover, on the experimental side, the signal-modulation problems first 
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discussed in Refs. [l6|, ll7| remain, to our knowledge, unsolved, while being of extreme 



importance on studying the feasibility of the experiment. 



Last, but not least, our findings should be compared with those in Ref. [19|, where the 
authors consider the cosmological evolution in a recently suggested new model of quantum 
initial conditions for the Universe. They find that the effective Friedmann equation incor- 
porates the effect of the conformal anomaly of quantum fields, and shows that their vacuum 
Casimir energy is completely screened and does not gravitate. 
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